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ABSTRACT 


Theoretical  estimates  of  high  frequency  sound  attenuation  and  backscattering  strengths 
from  suspended  sand  particles  in  the  ocean  and  in  an  estuary  are  given.  The  estimates 
are  obtained  using  the  measured  volume  concentrations  of  suspended  particle  in  the 
upper  layers  of  the  ocean  [Duykers,  1967],  from  the  particle  size  spectra  near  the  clear 
water  minimum  [McCave,  1984],  and  from  measurements  of  suspended  sediments  at 
different  depths  from  the  bottom  of  an  estuary,  [Thorne,  et  al,  1993].  Estimates  are  given 
for  frequencies  in  the  1-10  MHz  region.  Contributions  to  acoustic  attenuation  from 
scattering,  viscous  drag  and  thermal  conduction  are  evaluated  separately  using  existing 
theoretical  models.  Contribution  to  attenuation  from  thermal  conduction  is  negligible  at 
the  frequencies  used.  For  particles  in  the  sub  micrometer  to  10  \\m  range  attenuation  is 
predominantly  due  to  viscous  drag.  For  large  particles  scattering  is  the  main  cause  of 
attenuation.  Estimates  reveal  attenuation  and  backscattering  near  the  bottom  of  the 
estuary  is  quite  high.  At  5  MHz,  attenuation  at  a  depth  of  80  cm  and  10  cm  from  the 
bottom  are  approximately  0.1  dB/m  and  10  dB/m  respectively.  At  this  frequency 
backscattering  strength  at  80  cm  and  10  cm  depths  from  the  bottom,  are  approximately  - 
25  dB  and  -5  dB  respectively.  __ 
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Theoretical  Estimates  of  High  Frequency 
Acoustic  Attenuation  and  Backscattering  from 
Suspended  Sand  Particles  in  the  Ocean 
and  in  an  Estuary 


EXECUTIVE  SUMMARY 


This  study  looks  at  attenuation  and  backscattering  of  MHz  sound  by  suspended 
sand  particles.  It  is  part  of  a  project  to  investigate  the  use  of  sound  in  the  MHz 
region,  to  image  objects  in  sea.  The  study  examines  the  relative  importance  of 
various  mechanisms  that  causes  attenuation  of  sound  by  suspended  sand 
particles  and  gives  theoretical  estimates  of  attenuation  and  backscattering  at 
MHz  frequencies.  Estimates  are  given  from  regions  ranging  from  clear  water  to 
those  close  to  the  bottom  of  an  estuary.  The  study  shows  that  for  particles  that 
have  a  radius  less  than  10  microns  viscous  drag  is  the  predominant  mechanism 
for  attenuation,  while  for  large  particles  it  is  scattering.  Estimates  indicate  that 
attenuation  and  backscattering  close  to  the  bottom  of  the  sea  will  be  high,  which 
can  cause  problems  for  acoustic  visibility. 
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a  thermometric  conductivity 

Cp  specfic  heat  at  constant  pressure 

/ta  intensity  of  backsc  attered  sound  lm  from  a  unit  volume  of  scatterers 

/;  incident  intensity  of  sound 

BS  scattering  strength 


DSTO-TR-0078 


1.  Introduction 


The  importance  of  high  frequency  sound  attenuation  and  backscattering  from 
suspended  particulate  material  in  the  ocean  is  assessed  in  this  report,  using  theoretical 
models.  The  study  is  undertaken  to  evaluate  the  feasibility  of  using  high  frequency 
sound  in  the  MHz  region  to  obtain  an  image  of  an  object  in  sea.  One  of  the  problems 
associated  with  the  use  of  visible  light  is  the  large  scattering  caused  by  suspended 
particles. 

Most  of  the  early  work  in  the  distribution  of  total  suspended  material  (TSM)  came 
from  light  scattering  measurements.  Jerlov  [1953]  employing  light  scattering 
measurements  undertook  one  of  the  first  major  surveys  of  the  oceanic  TSM.  From  the 
data  he  was  able  to  identify  a  number  of  overall  trends  in  the  distribution  of  TSM, 
including  the  recognition  of  what  were  subsequently  termed  nepheloid  layers,  that  is 
layers  of  relatively  turbid  waters  that  can  extend  hundreds  of  metres  above  the  sea 
bed.  This  was  followed  by  a  series  of  investigations  by  Russian  scientists  in  the  1950's, 
[Lisitzin,  1972],  and  by  American  groups  in  the  1960's,  [Ewing  and  Thorndike,  1965; 
Connary  and  Ewing,  1972,  pi  69-1 84].  Recent  studies  of  the  global  distribution  of  the 
oceanic  TSM  have  been  carried  as  part  of  the  GEOSECS  and  HEBBLE  programs,  using 
sediment  traps  in  many  locations.  From  these  data  the  distribution  of  oceanic  TSM  can 
be  described  in  terms  of  a  three  layer  model  in  which  the  main  features  are,  a  surface 
layer,  a  clear  water  minimum  layer  and  a  deep  water  layer.  A  detailed  description  of 
the  distribution  of  TSM  is  given  in  Chester  [1990,  p321-345]. 

Attenuation  of  sound  by  suspended  material  can  be  due  to  scattering,  viscous  drag 
or  thermal  conduction.  Sections  2-4  describe  theoretical  evaluations  of  these  three 
mechanisms,  respectively,  Section  5  describes  the  method  to  evaluate  back  scattering 
strength.  In  section  6  numerical  estimates  of  attenuation  and  backscattering  strength 
are  given  for  4  cases.  The  first  case  uses  the  volume  concentration  of  particles  inferred 
for  the  upper  layers  (0-200  m)  of  the  ocean  by  Duykers  [1967],  using  low  frequency 
acoustic  data.  The  second  case  uses  the  particle  size  distribution  of  McCave  [1984], 
which  was  measured  near  the  clear  water  minimum.  The  third  case  involved  the  use 
of  the  measured  particle  size  distribution  in  the  Bay  of  Bengal,  in  the  Indian  ocean 
[Shifrin,1988,  p266].  Finally  estimates  are  obtained  from  suspended  sediment 
measurements  in  an  estuary,  at  different  distances  from  the  bottom  [Thome,  et  al, 
1993].  In  the  evaluation  of  estimates  particulate  matter  is  assumed  to  be  sand,  that  is 
silica.  Suspended  organic  matter,  such  as  phytoplanktons,  and  zooplanktons  are  poor 
scatterers  and  absorbers  when  compared  to  sand.  Estimates  are  given  for  frequencies 
in  the  1-10  MHz  region,  because  of  a  current  interest  in  mine  imaging,  for 
classification. 
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2.  Attenuation  Due  to  Scattering 


Preliminary  studies  on  the  use  of  available  models  to  estimate  scattering  of  sound  due 
to  sand  particles,  for  high  frequency  sound,  gave  the  following  results:  Application  of 
the  Rayleigh  approximation,  which  is  applicable  when  the  particle  size  is  much  less 
than  the  wavelength  of  sound,  for  particles  having  radius  larger  than  44  \xm  at 
frequencies  above  3  MHz  can  lead  to  significant  errors.  Results  from  the  non-rigid 
sphere  model  [Morse  and  Ingard,  1968,  p337],  are  always  less  than  the  rigid  sphere 
model  [Morse  and  Ingard,  1968,  p419]  by  a  factor  of  2  or  3,  while  the  inclusion  of  shear 
waves  has  effects  for  large  particles  (radius>150  \im).  Due  to  the  presence  of  such 
large  particles  near  the  sea  bottom  and  in  a  turbid  medium,  the  most  suitable 
approach  would  be  to  consider  the  suspended  material  to  be  non-rigid  spheres  and  to 
include  shear  waves. 


The  angular  scattered  intensity  function  |/(o,j)|  /  which  gives  the  angular 
distribution  of  the  scattered  intensity,  for  a  non  rigid  sphere  with  the  inclusion  of 
shear  waves  is  given  by  [Faran,  1951;  Hickling,  1962], 

| f(o,  Of  =  A  X  (2m  + 1)(2/  + 1) sin  r\m  sin  rj,  cos(rjm  -  T|,  )/>  (cos  0)Pm  (cos 6) 


(1) 


/,m= 0 


where  0  is  the  angle  between  o  and  i  ,  the  scattered  and  incident  directions  and  P,,Pm 
the  Legendre  polynomials  of  order  /  and  m  respectively,  k  is  the  wave  number  of  the 
2%f 

incident  wave(= - ),  f  the  frequency,  and  c  is  the  sound  speed  in  the  surrounding 

c 

medium.  Other  parameters  necessary  to  evaluate  the  angular  scattered  intensities  are 
given  below. 

=  tan  S„  (gKtan  <j)m  +  tan  qM  (s)) 
tan<t>m  +  tanPm(x) 

tmbn(x)  =  ^± 

yJx) 

tan  am  (x)  =  X^m+l  ^  -  m 

JJx) 

tanp  n(x)=W*^-m 

yjx) 

tanCmUPx2)  =  -0.  5x22  (~  Zl ) 


z,  = 


tan  am  ) 
tanam(x!)  +  l 


^2  = 


m  +m 


z,  = 


m  +m-l-0.5jc2  +tan  am(x2) 
_  m2  +m-0.5xl  +2tanam(^) 
tanam  (*,)  +  ! 
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z4  =  (tan<xm(x2)+l)z2 

tan(|)m  =^tanCm(JC„x2) 

Pi 

where  x  =  ka,  xx  =  bfl,  x2  =  k2a. 

O'trf  Oirf 

Wave  numbers  A:,  and  k2  are  — —  and  — — ,  respectively,  where  cx,  c2  are  the 

cx  c2 

velocity  of  compression  and  shear  sound  waves  in  the  particle  of  radius  a.  p ,  p  j  are 
the  densities  of  the  medium  and  particle  respectively.  jm  are  spherical  Bessel  functions 
of  the  first  kind  of  order  m  and  ym  are  spherical  Bessel  functions  of  the  second  kind  of 
order  m.  The  functions  ym  are  also  referred  to  as  spherical  Neumann  functions  of 
order  m,  [Morse  and  Ingard,  1968,  p33 7]. 

The  total  scattering  cross  section  of  a  particle  is  given  by: 

2rc[|/(d,0|2  sin  0d0 
o 

This  is  the  total  scattered  radiation  by  a  single  particle  per  unit  volume.  For  a  non  rigid 
sphere  with  the  inclusion  of  shear  waves,  using  eqn(l)  gives: 

X,=T?S(2m  +  1)sin2Tlm  (2) 

£  m=0 

Suppose  the  particle  number  density  (the  number  of  particles  per  m3)  is  N,  then  the 
total  scattering  loss  in  dB/m  is  given  by 

4.342AS,  (3) 

The  evaluation  of  the  spherical  Bessel  functions  jm,  of  all  orders  is  carried  out  using 
the  procedure  due  to  Miller,  [Abramowitz  and  Stegun,  1970,  p452].  The  spherical 
Bessel  functions  ym  are  obtained  by  first  evaluating  jn  for  negative  m  using 

ym_i  Qt) + jm+l  (x)  =  (2m  +  l)x~1jm  (x)  recursively,  and  then  calculating  ym  from  the 
relationship  yn(x)  =  (-l)m+1  j_m_x  (jc),  [Abramowitz  and  Stegun,  1970,  p439]. 


3.  Viscous  Attenuation 


This  section  examines  the  effect  due  to  viscosity  on  sound  attenuation.  Lamb  [1945] 
derived  the  theory  to  obtain  attenuation  of  sound  energy  by  particles  due  to  the 
viscosity  of  the  medium.  The  particles  were  considered  to  be  rigid  and  incompressible 
but  free  to  move.  This  was  tested  by  Urick  [1948]  who  obtained  good  agreement  with 
the  theory  of  Lamb  in  his  measurements  on  sand  and  kaolin  particles  at  megacycle 
frequencies.  Later  Hovem  [1980]  also  found  good  agreement  in  suspensions  provided 
the  volume  concentration  of  the  suspended  material  is  no  more  than  a  few  percent.  In 
the  ocean  volume  concentrations  of  suspended  particles  are  less  than  a  percent.  The 
theory  of  Lamb  is  therefore  suitable  to  study  viscous  attenuation  in  the  ocean. 
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Lamb  [1945]  calculated  the  attenuation  due  to  viscosity  by  finding  the  velocity 
potential  of  the  waves  scattered  by  the  particle  and  then  obtaining  the  average  rate  at 
which  work  is  being  done  over  the  spherical  surface  surrounding  the  particle.  For 
ka<\,  k  being  the  wave  number  of  incident  acoustic  wave  ,  he  showed  that 
attenuation  due  to  viscosity  becomes  important  when  the  radius  of  the  particle  is  less 

,  2v  - 

than  or  comparable  to  P  —  ( — )2,  where  V  is  the  kinematic  viscosity  of  the  fluid  and 

C0(=  2%f)  is  the  radian  frequency.  For  a  «  p_1  or  a  »  p_1  viscous  effects  are  minimal, 
[Urick  and  Ament,  1949].  The  attenuation  due  to  viscosity  in  dB/m  for  particles  of 
radius  a  is  given  by,  [Lamb,  1945,  p660] 

4:c 

4.342 A— Real  (iA.B,  -  iA0B0 ) 


where  N  is  the  number  of  particles  per  m3 , 
A  -(3+3iha-(ha)2)(ka)3eii" 

1  (ka)2(l  +  iha)  +  (2  +  2ika-(ka)2)(ha)2 

B  Pi-3pp wika) 

Pi-Po-PoW/iC^Mi 


y(ka)  = 


sin  (ka) 
(ka)3 


co  s(ka) 
(ka)2 


(4) 


f1(ka)  =  e~ika(- 


-3  i 
(ka)3 


3 

(ka)4  + 


(ka)2) 


ha  =  (l-i)$a. 

Aq  and  B0  are  the  values  of  Ax  and  B1  when  ha  -» 
limited  to  extremely  small  values  of  ka. 


00 .  The  above  expression  is  not 


A  much  simpler  formula  used  by  Urick  [1948]  and  also  by  Ahuja[1970]  for  rigid 
particles  of  radius  'a'  and  valid  for  ka  « 1  is, 

4.342  Ar(-7M3)*(£i--l)- - - - 

3  P  s2  +  (^  +  f)2 
P 


9  1 

5  4p<3  1  +  pa 


t  =  0.5  + 


9 

4p  a 
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4.  Thermal  Attenuation 


In  this  section  the  effect  of  thermal  conduction  on  sound  attenuation  is  examined. 
Epstein  and  Carhart  [1953],  in  their  investigation  of  sound  absorption  by  fog  particles, 
obtained  a  result  to  estimate  attenuation  of  sound  due  to  thermal  conduction.  Allegra 
and  Hawley  [1971]  found  good  agreement  with  their  attenuation  results  in  emulsions 

AkFg 

where  thermal  transport  processes  are  important.  Provided  — —  «1,  where  O  is  the 

c 

X 

thermometric  conductivity  of  the  medium,  defined  by  - ,  T  being  the  thermal 

pcp 

conductivity  and  Cp  the  specific  heat  at  constant  pressure,  an  estimate  of  attenuation 

due  to  thermal  conduction  can  be  obtained.  For  frequencies  of  interest  here  (1-10 
MHz),  and  for  sea  water  this  ratio  is  «1.  The  thermal  attenuation  in  dB/m  due 
particles  of  radius  a  is  then  given  by  Epstein  and  Carhart  [1953]  as 

4.342(4tcW(y  - 1)(— )a(l -If  Real  (— ) ),  (5) 

c  F 


where 


ao 

,  „  ,  1  *  sink 

and  F  —  ( - )  —  y - ! - 

1  -ib  sinbj-bjcosfy 

2  nf  - 

The  parameter  b  is  defined  by,  b  =  (l+0( — —)2a.  The  suffix  1  indicates  particle. 

2a 

x 

X  =  — ,  and  a  is  the  coefficient  of  volume  expansion.  The  expression  is  valid  for 
ka  « 1.  For  other  values  of  ka  attenuation  from  thermal  conduction  is  negligible.The 


necessary  thermal  data  for  sand  and  sea  water  are  obtained  from  Kaye  and  Laby 
[1972],  and  Walton  Smith  [1974]. 


5.  Backscattering  Strength 

Backscattering  strength (BS)  in  dB  is  defined  as, 

BS  =  10Logl„(!f-), 

■** 

where  /fe  is  the  intensity  of  the  backscattered  sound  at  a  distance  of  lm  from  a  unit 
volume  of  scatterers  and  /,.  is  the  incident  intensity.  If  scattering  is  omnidirectional 
then 
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vrv 

BS  =  lOLog(-^-).  (6) 

4k 

This  would  be  the  case  if  ka  « 1.  Otherwise 

BS  =  1OL06„(W|/(18O-)|2)  (7) 

where  |/(180‘)|  is  the  angular  scattered  intensity  function  given  by  eqn(l),  evaluated 

for  0  =  180°.  The  angular  scattered  intensity  function  is  the  same  as  the  commonly 
described  volume  scattering  function. 


6.  Numerical  Estimates 


In  this  section  numerical  estimates  for  attenuation,  and  backscattering  strength  due  to 
suspended  sand  particles  are  presented  from  experimentally  measured  particle 
concentrations  or  distributions.  Figurel(a)  shows  the  total  attenuation  and 
contributions  due  to  scattering,  viscous  drag  and  thermal  conduction,  for  seawater 
containing  particles  with  a  uniform  radius,  as  a  function  of  particle  radius  at  5  MHz. 
The  volume  concentration  of  particles  is  taken  as  6x10”*,  this  being  inferred  by 
Duykers  [1967]  from  acoustic  measurements  in  the  upper  layers  (0  -200m)  of  the 
ocean.  The  particle  radius  was  varied  from  0.5  pm  to  50  | lm.  Attenuation  due  to 
thermal  conduction  is  negligible  at  all  particle  sizes.  For  particles  with  radii  less  than 
10  pm  viscous  drag  is  the  dominant  effect  of  attenuation,  while  beyond  10  [lm 
scattering  contributes  the  most  to  attenuation,  [Uric,1948].  The  variation  of  viscous 
attenuation  with  radius  is  explained  as  follows. 

When  sound  wave  passes  through  the  medium,  the  particles  will  oscillate  with  the 
fluid.  If  the  medium  is  viscous  then  there  will  be  loss  of  energy  due  to  drag.  The  drag 
force  that  causes  dissipation  depends  on  the  relative  speed  between  particle  and  fluid 
and  on  the  surface  area  of  the  particle.  Small  particles  have  a  low  inertia  and  vibrate  to 
and  fro  along  with  the  fluid.  Hence  the  relative  speed  is  small.  On  the  other  hand  large 
particles  have  a  high  inertia  and  tend  to  lag  behind.  However  for  a  constant  volume 
concentration  of  particles,  the  available  surface  area  will  decrease  with  increase  in 
radius.  The  combination  of  these  two  effects  results  in  viscous  attenuation  at  a 
particular  frequency  being  important  only  for  smaller  radii.  This  is  seen  in  Fig.l(a). 
Fig.  1(b)  shows  the  backscattering  strength  at  5  MHz  as  a  function  of  the  particle 
radius.  The  backscattering  increases  with  particle  radius. 

From  measurements  of  particle  size  spectra  near  the  clear  water  minimum,  McCave 
[1984]  showed  that  the  particle  number  data  can  nearly  always  be  fitted  by  a  power 
law  distribution  over  a  large  part  of  the  measured  range.  Expressed  as  a  cumulative 
number  N0  as  a  function  of  particle  diameter  d,  the  size  distribution  is  described  by 
N  0=  kcd~3 

for  d  in  the  range  1-100  pm,  and  kc  =  5X109per  m3.  N0  is  the  number  of  particles 
with  diameter  >  d .  This  indicates  equal  particle  volumes  in  logarithmically  increasing 
size  grades.  Fig.  2(a)  gives  the  attenuation  due  to  viscous  drag,  scattering  and  the  total, 
for  frequencies  1-10  MHz  for  such  a  distribution.  The  attenuation  values  are  small. 
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Fig.2(b)  shows  the  backscattering  strength  for  frequencies  1-10  MHz  for  such  a 
distribution.  Backscattering  strength  increases  with  frequency. 

Shiffin  [1988,  p266]  has  presented  data  on  the  distribution  of  particles  obtained 
from  measurements  made  in  the  summer  of  1979  in  the  Indian  Ocean.  One  of  them  is 
the  result  from  the  Bay  of  Bengal  in  the  vicinity  of  the  Ganges.  The  particle 
distribution  given  is  obtained  at  depth  of  10m.  Most  of  the  particles  have  radii  less 
than  1  \im .  Figure  3(a)  gives  the  attenuation  due  to  viscous  drag,  scattering  and  total, 
for  frequencies  1-10  MHz.  Attenuation  due  to  thermal  conduction  is  negligible.  Most 
of  the  attenuation  is  due  to  viscous  drag,  which  arises  due  to  the  presence  of  a  large 
number  of  particles  that  have  radii  less  than  1  |im.  Figure  3(b)  shows  the 
backscattering  strength  for  frequencies  1-10  MHz.  It  increases  with  frequency  but  the 
values  are  similar  to  that  of  clear  water  due  to  the  absence  of  large  particles  that 
contribute  to  scattering. 

Thome  et  al  [1993],  carried  out  suspended  sediment  measurements  in  an  estuary  at 
different  depths  close  to  the  bottom.  The  bed  consists  of  well-sorted  fine  sand  with 
little  clay  or  silt  present  in  the  bed  or  in  the  water  column.  Analysis  of  the  samples  at 
10  and  80  cm  above  the  bed  showed  a  Gaussian  probability  distribution  with  an 
average  radius  of  85  \im  and  78  \un  respectively  and  a  standard  deviation  of  20  \im. 
The  total  concentrations  of  particles  at  these  two  depths  are  1  kg  /  m3  and  10  2  kg  /m3. 
Fig.  4(a)  shows  the  attenuation  due  to  viscous  drag  and  scattering  80  cm  from  the 
bottom  for  frequencies  in  the  1-10  MHz  region,  while  Fig.  4(b)  is  for  measurements  10 
cm  from  the  bottom.  Due  to  the  large  sediment  particle  size  present  viscous 
attenuation  is  small  compared  to  scattering.  Figures  5(a)  and  5(b)  are  the 
backscattering  strength  as  a  function  of  frequency  for  the  2  depths  80  cm  and  10  cm 
from  the  bottom.  The  plateau  observed  in  the  backscattered  strength  is  due  to  the 
modes  of  vibration  of  the  elastic  sphere  being  excited  at  certain  frequencies  which  then 
interfere  destructively  with  the  backscattered  signal,  [Thome  et  al  ,1992].  It  is  a 
resonance  phenomenon.  It  is  absent  if  the  particles  are  considered  rigid.  We  are  able  to 
observe  this  in  the  estuarine  measurements  due  to  the  presence  of  large  particles,  but 
it  is  absent  in  the  region  of  the  clear  water  minimum. 


7 


8ACKSCATTERING  STRENGTH(dB)  LOG(ATTENUATION  IN  dB/m) 


BACKSCATTERING  STRENGTH(dB  ) 


BACKSCATTERING  STRENGTH(dB) 


DSTO-TR-0078 


2.5  5.0  7.5  10.0 


FREQUENCY  (MHz) 


2.5  5.0  7.5  10.0 


FREQUENCY  (MHz) 


10 


Figure  3:  (a)  Total  attenuation  and  attenuation  due  to  scattering,  viscous  drag  and  thermal 
conduction,  in  dB/m  and  (b)  Backscattering  strength,  in  the  Bay  of  Bengal  in  the  Indian  ocean, 
as  a  function  of  frequency. 
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7.  Conclusion 


This  paper  provides  theoretical  estimates  of  attenuation  and  backscattering  strengths 
due  to  suspended  sand  particles  at  frequencies  in  the  1-10  MHz  region.  Four  different 
situations  are  examined.  The  first  case  involved  the  use  of  volume  concentration  of 
particles  in  the  shallow  layers  (0-200  m)  of  the  ocean.  Results  showed  that  thermal 
attenuation  is  negligible  for  all  particle  sizes  (<1(T3  dB/m  at  5  MHz).  For  particles 
with  radii  less  than  10  pm  viscous  drag  is  the  dominant  cause  of  attenuation.  For 
particles  having  radii  near  1  | lm  viscous  attenuation  is  approximately  10-1  dB/m  at  5 
MHz.  It  decreases  with  radius  beyond  0.5  pm.  Attenuation  due  to  scattering  for 
particles  having  radii  near  1  pm  is  less  than  10-5  dB/m  at  5  MHz.  Attenuation  due  to 
scattering  increases  with  radius  and  beyond  10  p m  it  is  the  main  contribution  to 
attenuation.  For  particles  near  50  pm  attenuation  due  to  scattering  is  10_1  dB/m  at  5 
MHz  while  viscous  attenuation  is  about  10-3  dB/m.  Backscattering  strength  increases 
with  particle  radius  and  varies  from  -80  dB  to  -25  dB  as  particle  radius  changes  from 
0.5  to  50  pm. 

The  second  case  involved  the  use  of  the  particle  size  distribution  of  McCave  [1984], 
valid  near  the  clear  water  minimum.  This  distribution  is  obtained  from  a  large  number 
of  measurements  in  the  North  West  Atlantic,  Mediterranean  and  the  East  Pacific.  For 
this  distribution,  in  the  1-5  MHz  region  the  total  attenuation  is  mainly  due  to  viscous 
drag  and  has  a  value  of  approximately  10^*  dB/m.  Beyond  5  MHz  scattering  too 
becomes  important  and  the  total  attenuation  is  5X1 dB/m.  Viscous  attenuation  is 
fairly  constant  over  1-10  MHz  frequency  range,  while  attenuation  due  to  scattering 
varies  from  approximately  1(T7  to  10*4  dB/m  over  this  frequency  range. 
Backscattering  strength  increases  with  frequency  and  varies  from  nearly  -80  to  -50  dB 
over  the  frequency  range  1  to  10  MHz. 

The  third  case  involved  the  measured  particle  size  distribution  near  the  Bay  of 
Bengal  in  the  Indian  Ocean,  at  10  m  depth.  The  radii  of  all  particles  were  less  than 
10  pm .  The  attenuation  is  mainly  due  to  viscous  drag  and  has  a  near  constant  value 

of  approximately  4xl0-3  dB/m  over  the  1-10  MHz  frequency  range.  Backscattering 
strength  increases  with  frequency  and  shows  a  similar  variation  to  that  obtained  for 
clear  water. 

The  last  case  used  suspended  sediment  measurements  made  in  an  estuary  at  two 
different  levels  80  cm  and  10  cm  from  the  bottom.  At  both  depths  the  main 
contribution  of  attenuation  is  from  scattering,  with  viscous  attenuation  being 
negligible  in  comparison.  This  arises  due  to  the  presence  of  large  particles  present  at 
these  depths,  the  average  radius  being  around  80  pm.  For  particles  with  such  radii 
viscous  drag  is  small.  At  80  cm  from  the  bottom,  attenuation  varies  from  10"3  dB/m  to 
about  10"1  dB/m  over  the  frequency  range  1  to  5  MHz.  In  the  5-10  MHz  region 
variation  is  small,  attenuation  being  close  to  10-1  dB/m.  At  10  cm  from  the  bottom. 
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attenuation  varies  from  10-1  dB/m  to  10  dB/m  over  the  frequency  range  1  to  5  MHz, 
while  in  the  5-10  MHz  region  it  is  nearly  constant  around  10  dB/m.  Backscattering 
strengths  at  80  cm  from  the  bottom  varies  from  -42  dB  to  -25  dB  over  the  frequency 
range  1-5  MHz,  remaining  fairly  constant  at  close  to  -25  dB,  in  the  5-10  MHz  range.  At 
10  cm  above  the  bottom,  backscattering  strength  varies  from  -22  dB  to  -5  dB  in  the  1-5 
MHz  region,  remaining  constant  at  around  -5  dB  in  the  5-10  MHz  range. 

The  study  shows  that  the  predominant  mechanism  for  attenuation  is  viscous  drag 
for  particles  with  radii  less  than  10 1 lm,  while  those  with  larger  radii  it  is  scattering. 
Although  attenuation  results  in  loss  of  acoustic  energy,  backscattering  is  of  more 
serious  concern.  It  can  cause  glitter  in  the  screen  and  can  mask  the  image  of  the  object 
that  is  of  interest.  Large  sand  particles  are  mostly  present  close  to  the  sea  bottom  and  it 
is  here  that  acoustic  visibility  would  pose  problems. 


14 


DSTO-TR-0078 


8.  Acknowledgments 


The  author  wishes  to  thank  Dr.  I.  S.  F.  Jones  and  Dr.  P.  Mulheam  for  their  comments 
and  suggestions. 


References 


1.  Abramowitz,  M.,  and  I.  A.  Stegun,  Handbook  of  Mathematical  functions,  Dover 
publications.  New  York,  (1970). 

2.  Ahuja,  A.  A.,  Effect  of  particle  viscosity  on  propagation  of  sound  in  suspensions  and 
emulsions,  J.  Acous.  Soc.  Am.,  51, 182-191  (1972). 

3.  Allegra,  J.  R.,  and  S.  A.  Hawley,  Attenuation  of  sound  in  suspensions  and  emulsions: 
Theory  and  Experiments,  J.  Acous.  Soc.  Am.,  51, 1545-1564  (1972). 

4.  Chester,  R.,  Marine  Geochemistry,  Uniwan  Hyman  Ltd,  London,  (1990). 

5.  Connary,  S.  C.,  and  M.  Ewing,  The  nepheloid  layer  and  bottom  circulation  in  the 
Guinea  and  Angola  Basins,  Studies  in  Physical  Oceanography,  A.L.  Gordon  (ed),  169- 
84,  London:  Gordon  and  Breach,  (1972). 

6.  Duykers,  L.  R.  B.,  Sound  attenuation  in  liquid-solid  mixtures,  J.  Acous.  Soc.  Am.,  41, 
1330-1335  (1967). 

7.  Epstein,  S.  P.,  and  R.  R.  Carhart,  The  absorption  of  sound  in  suspensions  and 
emulsions.  1.  Water  fog  in  air,  J.  Acous.  Soc.  Am.,  25, 553-565  (1953). 

8.  Ewing,  M.,  and  RThomdike,  Suspended  matter  in  deep  ocean  water,  Science,  147, 
1291-4  (1965). 

9.  Faran,  J.  J.,  Sound  Scattering  by  solid  cylinders  and  spheres,  J.  Acous.  Soc.  Am.,  23, 
405-418  (1951). 

10.  Hickling,  R.,  Anaylsis  of  echoes  from  a  solid  elastic  sphere  in  water,  J.Acoust.  Soc.  Am., 
34, 1582-1592  (1962). 

11.  Hovem,  Jens  M.,  Viscous  attenuation  of  sound  in  suspensions  and  high  porosity  marine 
sediments,].  Acous.  Soc.  Am.,  67, 1559-1563  (1980). 


15 


DSTO-TR-0078 


12.  Jerlov,  N.  G.,  Particle  distribution  in  the  ocean,  Rep.  Swed.  Deep-Sea  Exped.,  3,  71- 
97(1953). 

13.  Kaye,  G.  W.  C.,  and  T.  H.  Laby,  Tables  of  physical  and  chemical  constants  and  some 
mathematical  functions,  Longman  Press,  New  York,  (1972). 

14.  Lamb,  H.,  Hydrodynamics,  Dover  Publications,  New  York,  361-363  (1945). 

15.  Lisitzin,  A.  P.,  Sedimentation  in  the  world  ocean,  Soc.  Econ.  Paleo.  Min.  Spec.Publ., 
no  17,  (1972). 

16.  McCave,  I.  N.,  Size  spectra  and  aggregation  of  suspended  particles  in  the  deep  ocean, 
Deep  Sea  Res.,  31, 329-352  (1984). 

17.  Morse,  P.  M.,  and  K  U.  Ingard,  Theoretical  Acoustics,  McGraw-Hill  Inc,  New  York 
(1968). 

18  Shifrin,  K.  S.,  Physical  Optics  of  Ocean  Water,  Amercan  Institute  of  Physics,  New 
York  (1988). 

19.  Thome,  P.  D.,  L.  Hayhurst,  and  V.  L.  Humphrey,  Scattering  by  nonmetalic  spheres, 
Ultrasonics.,30, 15-20  (1992). 

20.  Thome,  P.  D.,  P.  J.  Hardcastle,  and  R.  L.  Soulsby,  Analysis  of  acoustic  measurements 
of  suspended  sediments,  J.  Geophy.  Res.,  98, 899-910  (1993). 

21.  Uric,  R.  J.,  The  absorption  of  sound  in  suspensions  of  irregular  particles,  J.  Acous.  Soc. 
Am.,  20, 283-289  (1948). 

22.  Urick,  R.  J.,  and  W.  S.  Ament,  The  propagation  of  sound  in  composite  media,  J.  Acous. 
Soc.  Am.,  21, 115-119  (1949). 

23.  Walton  Smith,  F.  G.,  Handbook  of  Marine  Sciences,  CRC  press,  (1974). 


16 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


UNCLASSIFIED 


REPORT  NO.  ARNO.  REPORT  SECURITY  CLASSIFICATION 

DSTO-TR-0078  AR-008-946  UNCLASSIFIED 


TITLE 

Theoretical  Estimates  of  High  Frequency  Acoustic  Attenuation  and  Backscattering  from  Suspended  Sand 
Particles  in  the  Ocean  and  in  an  Estuary 


AUTHOR(S) 

R.A.  Thuraisingham 

CORPORATE  AUTHOR 

Aeronautical  and  Maritime  Research 

Laboratory 

GPO  Box  4331 

Melbourne  Victoria  3001 

REPORT  DATE 

TASK  NO. 

SPONSOR 

November  1994 

NAV  92/257 

Navy 

FILE  NO. 

REFERENCES 

PAGES 

510/207/0176 

23 

26 

CLASSIFICATION/LIMITATION  REVIEW  DATE 

CLASSIFICATION/RELEASE  AUTHORITY 

Maritime  Operations  Division 

SECONDARY  DISTRIBUTION 


Approved  for  public  release 


ANNOUNCEMENT 


Announcement  of  this  report  is  unlimited 


KEYWORDS 

High  frequencies  Attenuation  particle  size  Acoustic  attenuation 

Sound  scattering  Viscous  drag  Backscattering 


ABSTRACT 

Theoretical  estimates  of  high  frequency  sound  attenuation  and  backscattering  strengths  from  suspended  sand 
particles  in  the  ocean  and  in  an  estuary  are  given.  The  estimates  are  obtained  using  the  measured  volume 
concentrations  of  suspended  particle  in  the  upper  layers  of  the  ocean  [Duykers,  1967],  from  the  particle  size 
spectra  near  the  clear  water  minimum  [McCave,  1984],  and  from  measurements  of  suspended  sediments  at 
different  depths  from  the  bottom  of  an  estuary,  [Thome,  et  al,  1993].  Estimates  are  given  for  frequencies  in  the 
1-10  MHz  region.  Contributions  to  acoustic  attenuation  from  scattering,  viscous  drag  and  thermal  conduction 
are  evaluated  separately  using  existing  theoretical  models.  Contribution  to  attenuation  from  thermal  conduction 
is  negligible  at  the  frequencies  used.  For  particles  in  the  sub  micrometer  to  10  |lm  range  attenuation  is 
predominantly  due  to  viscous  drag.  For  large  particles  scattering  is  the  main  cause  of  attenuation.  Estimates 
reveal  attenuation  and  backscattering  near  the  bottom  of  the  estuary  is  quite  high.  At  5  MHz,  attenuation  at  a 
depth  of  80  cm  and  10  cm  from  the  bottom  are  approximately  0.1  dB/m  and  10  dB/m  respectively.  At  this 
frequency  backscattering  strength  at  80  cm  and  10  cm  depths  from  the  bottom,  are  approximately  -25  dB  and 
-5  dB  respectively. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 
UNCLASSIFIED 


Theoretical  Estimates  of  High  Frequency  Acoustic  Attenuation  and 
Backscattering  from  Suspended  Sand  Particles  in  the  Ocean  and  in  an  Estuary 

R.A.  Thuraisingham 

DSTO-TR-0078 

DISTRIBUTION  LIST 


Director,  AMRL 

Chief  Maritime  Operations  Division 

Author  Ranjit  A.  Thuraisingham  3  copies 

Library,  AMRL  Maribymong 
Library,  AMRL  Fishermens  Bend 

Chief  Defence  Scientist  (for  CDS,  FASSP,  ASSCM)  1  copy  only 

Director,  ESRL 

Head,  Information  Centre,  Defence  Intelligence  Organisation 
OIC  Technical  Reports  Centre,  Defence  Central  Libraiy 

Officer  in  Charge,  Document  Exchange  Centre  8  copies 

Senior  Defence  Scientific  Adviser 

Air  Force  Scientific  Adviser,  Russell  Offices 

Scientific  Adviser  -  Policy  and  Command 

Senior  Librarian,  Main  Library  DSTOS 

Librarian,  DSD,  Kingston  ACT 

Serials  Section  (M  List),  Deakin  University  Library,  Deakin  University,  Geelong  3217 
NAPOC  QWG  Engineer  NBCD  c/-  DENGRS-A,  HQ  Engineer  Centre,  Liverpool 
Military  Area,  NSW  2174 

ABCA,  Russell  Offices,  Canberra  ACT  2600  4  copies 

Librarian,  Australian  Defence  Force  Academy 

Head  of  Staff,  British  Defence  Research  and  Supply  Staff  (Australia) 

NASA  Senior  Scientific  Representative  in  Australia 
INSPEC:  Acquisitions  Section  Institution  of  Electrical  Engineers 
Head  Librarian,  Australian  Nuclear  Science  and  Technology  Organisation 
Senior  Librarian,  Hargrave  Library,  Monash  University 
Library  -  Exchange  Desk,  National  Institute  of  Standards  and  Technology,  US 
Acquisition  Unit  (DSC-EO/GO),  British  Library,  Boston  Spa,  Wetherby,  Yorkshire  LS23 
7BQ,  England 

Library,  Chemical  Abstracts  Reference  Service 
Engineering  Societies  Library,  US 

Documents  Librarian,  The  Center  for  Research  Libraries,  US 
Navy  Scientific  Adviser 

Army  Scientific  Adviser,  Russell  Offices  -  data  sheet  only 
Director  General  Force  Development  (Land)  -  data  sheet  only 
DASD,  APW2-1-OA2,  Anzac  Park  West,  Canberra  ACT  -  data  sheet  only 
SO  (Science),  HQ  1  Division,  Milpo,  Enoggera,  Qld  4057  -  data  sheet  only 
Librarian  -  AMRL  Sydney  -  data  sheet  only 

Counsellor,  Defence  Science,  Embassy  of  Australia  -  data  sheet  only 
Counsellor,  Defence  Science,  Australian  High  Commission  -  data  sheet  only 
Scientific  Adviser  to  DSTC  Malaysia,  c/-  Defence  Adviser  -  data  sheet  only 
Scientific  Adviser  to  MRDC  Thailand,  c/-  Defence  Attache  -  data  sheet  only 

Dr  Ian  Jones,  MOD  Sydney 
Mr  Jim  Thompson,  MOD  Sydney 
Dr  Phillip  Mulheam,  MOD  Sydney 
Dr  David  Blair,  MOD  Sydney 

Mr  Ashley  Deacon,  Business  Manager  -  Sensors,  GEC  Mavconi  Systems  Pty  Ltd 
Faraday  Park,  Railway  Road,  Meadowbank,  NSW  2114 
Dr  D.  Robinson,  Ultrasonics,  CSIRO  Division  of  Radiophysics,  126  Greville  Street, 

Chatswood,  NSW  2067 


